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Abstract 

We construct a low-energy effective field theory that permits the complete treatment of 
isospin-breaking effects in nonleptonic weak interactions to next-to-leading order. To this 
end, we enlarge the chiral Lagrangian describing strong and AS" = 1 weak interactions 
by including electromagnetic terms with the photon as additional dynamical degree of 
freedom. The complete and minimal list of local terms at next-to-leading order is given. 
We perform the one-loop renormalization at the level of the generating functional and 
specialize to K — * decays. 
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1 Introduction 



Although isospin violation in nonleptonic weak interactions has been investigated many 
times in the past systematic treatments have appeared only rather recently |1|, |2], ||, [| [| . 
The topic is both of general interest and of considerable phenomenological relevance. 
Precise determinations of weak decay amplitudes are needed for many purposes, in par- 
ticular for a reliable calculation of CP violation in the K° — K° system. In the standard 
model, isospin violation arises from the quark mass difference m u — and from electro- 
magnetic corrections. Although these effects are expected to be small in general they are 
amplified in nonleptonic weak transitions. Because of the suppression of amplitudes with 
AI > 1/2, isospin violation in the dominant AI = 1/2 amplitudes leads to significantly 
enhanced corrections for the sub-dominant amplitudes. In fact, a quantitative analysis of 
the AI = 1/2 rule is only possible with the inclusion of isospin- violating effects. 

At first order in a systematic low-energy expansion, isospin breaking in the leading 
octet amplitudes of nonleptonic kaon decays is of order G 8 (m u — m d ) and G$e 2 where G 8 
denotes the strength of the effective octet coupling. The corrections appear in the mass 
differences of charged and neutral mesons, via tc° — rj mixing and through electromagnetic 
penguins in the effective nonleptonic weak Hamiltonian. However, there are good 
reasons to believe that the problem cannot be understood at lowest order only. For 
instance, the resulting (tree- level) corrections do not produce a AI = 5/2 component for 
which there is some phenomenological evidence [p]]. 

The chiral realization of isospin violation due to the light quark mass difference is 
available also at next-to-leading order. The purpose of this paper is to close the gap in 
the electromagnetic sector by 

• completing the construction of the effective chiral Lagrangian of 0(G%e 2 p 2 ) and 

• performing the complete renormalization at the one-loop level for nonleptonic weak 
transitions including electromagnetic corrections. 

As our notation indicates, we only consider corrections to the leading octet part of the 
nonleptonic weak Hamiltonian. The results are applicable to the analysis of both K — > 2n 
[|, |, |, and K -> 3vr decays. 

We start in Sec. |2| by recalling the ingredients for the construction of effective theories 
of strong, electromagnetic and nonleptonic weak interactions. In Sec. ^ we review the 
effective Lagrangian of lowest order. For this Lagrangian, we evaluate the one-loop diver- 
gence functional by standard heat-kernel techniques in Sec. |j. The new parts are terms of 
0(Gse 2 p 2 ) which arise also from using the equations of motion to transform to the stan- 
dard bases for the nonleptonic weak Lagrangian of 0{G^,p i ) || and for the electromagnetic 
Lagrangian of 0(e 2 p 2 ) 0. In the following section we construct the complete and minimal 
Lagrangian of 0(G$e 2 p 2 ) making use of CPS symmetry [|IU], Cayley-Hamilton relations, 
partial integration in the action and of the equations of motion. We order the terms in the 
effective Lagrangian according to their physical relevance: K — > tttt amplitudes receive 
contributions from the first 12 operators, the next two appear in K — > 3tt and the rest 
turns out not to be relevant phenomenologically. In Sec. ||, we present the divergences 
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for the three K — ► 2ir amplitudes and compare with the results of direct one-loop calcu- 
lations [/J. We summarize our findings in Sec. [7]. Various quantities appearing in the 
heat-kernel expansion of the generating functional are collected in the Appendix. 



2 Symmetries 

For a complete treatment of isospin-breaking effects in nonleptonic kaon decays, an ap- 
propriate effective Lagrangian with the pseudoscalar octet and the photon as dynamical 
degrees of freedom has to be used. The symmetries of the standard model are serving as 
the basic guiding principles for its construction. Starting with QCD in the chiral limit 
m u = m d — fn s — 0, the resulting symmetry under the chiral group G = SU (3) l X SU (3)r 
is spontaneously broken to SU(3)y- The pseudoscalar mesons (n,K,r]) are interpreted 
as the corresponding Goldstone fields <pi (i — 1, . . . , 8) acting as coordinates of the coset 
space SU(3)l x SU(3)r/ SU(3)y- The coset variables wl,,r(v?) are transforming as 

G 1 

Ur{<p) -> gRUR{<f)h(g,<f) , 
9 = (9l,9r) e SU{3) L xSU(3) R , (2.1) 

where h(g,(p) is the nonlinear realization of G 
The photon field is introduced in 



^Wr^ _ irp)u R - u\(d^ - il^)u L \ (2.2) 
by adding appropriate terms to the usual external vector and axial- vector sources v^, a^. 

In = v u ~ a u ~~ cQlA-u , 



r 



i-i 



Vft + an- zQrA^ . (2.3) 



The 3x3 matrices Ql,r are spurion fields with the transformation properties 

G t G i 

Ql -> 9lQl9l > Qfl -» 9rQr9r (2-4) 
under the chiral group. We also define 

Ql := u[Q L u L , Q fl := u R Q R u R (2.5) 

transforming as 

Q R $ h{g^)Q R h{g^Y l . (2.6) 
At the end, one identifies Ql,_r with the quark charge matrix 

(2.7) 





" 2/3 













-1/3 













-1/3 
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External scalar and pseudoscalar sources are combined in 



x = s + ip 



(2.8) 



For the construction of the effective Lagrangian, it is convenient to introduce the quan- 
tities 

X± = UrXu l ± u{x^u R (2.9) 

with the same chiral transformation properties as Ql, Qr in (|2.6|) . 

After integrating out the heavy degrees of freedom, the AS* = 1 weak interactions can 



be described in terms of an effective four-fermion Hamiltonian With respect to the 
chiral group G, this effective Hamiltonian transforms as the direct sum 



(8 l ,1r) + (27l,Ir) + (Sl,8 r ) 



(2.10) 



where the first piece, contributing only to AI = | transitions, is largely dominant. In 
this work we shall consider only the electromagnetic corrections induced by the dominant 
octet part of the effective Hamiltonian. To this end we introduce a weak spurion A that 
is finally taken at 

x . x T 

A 6 - zA 7 



A 
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(2.11) 



where Xqj are Gell-Mann matrices. In analogy to fl2.5|) we also define 

A := u{\u L , 



(2.12) 

transforming again as in ( |2.(j| ) under chiral transformations. 

Although CP is broken by the weak interactions, the AS* = 1 transitions are still 
invariant under the so-called CPS symmetry |l(J: a CP transformation followed by a 
subsequent interchange of d and s quarks. This symmetry is also obeyed by strong and 
electromagnetic interactions, provided the 2-3 indices of the external fields are also ex- 
changed appropriately (this implies, in particular, the exchange m s <-> in the mass 
terms). The explicit CPS transformation properties of the several building blocks intro- 
duced so far are given by 



with 
and 



u^(x) C -^> S —e(/j)Su^(x)S 

x ±(x) C ^ S ±Sxi(x)S , 

QlM °™ SQ T LtR (x)S, 

A(x) C 4 5 SA T (S)S , 



x = (x°, - x) , e(0) = 1 , e(l) = e(2) = e(3) = -1 



S 



1 
1 
1 



(2.13) 
(2.14) 

(2.15) 
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3 The effective Lagrangian at lowest order 



With the building blocks introduced in the previous section we may now assemble our 
effective Lagrangian. We adopt an expansion scheme where the n-th order is related to 
terms of order p n in the strong and weak sector and to terms of order e 2 p n ~ 2 in the 
electromagnetic sector where p denotes a typical meson momentum. Terms of 0(e 4 ) will 
be neglected throughout. 

To lowest order (n — 2), our effective theory consists of the following parts: the strong 
sector is represented by the nonlinear sigma model in the presence of the external sources 
v ni a n,X HH an d the photon coupling introduced in fl2T~ 

F 2 
1 



— («X + X+> • (3-1) 



The symbol ( ) denotes the trace in three-dimensional flavour space and F is the pion 
decay constant in the chiral limit. Explicit chiral symmetry breaking by the non- vanishing 
masses of the light quarks is achieved by evaluating the generating functional at 



X = 2BM quaA = 2B 



m u 
m d 
m, 



(3.2) 



The quantity B is related to the quark condensate in the chiral limit by (0|gg|0) = —F 2 B. 

The (8l, 1r) piece of the nonleptonic weak interactions is represented by the well- 
known Cronin Lagrangian fL4|. 



F 2 , 2 = G 8 F 2 A + h.c. . (3.3) 



At lowest order, the parameter G 8 can be determined [|HJ from K — > 2n decays to be 
\G 8 \ ~ 9 x 10- 6 GeV" 2 ~ 5(G F /V2)\V ud V us \. 

Now also the electromagnetic interaction has to be included. Apart from the necessary 
modification in (|2.2|) , we have to add a kinetic term for the photon field, 

- \f^ v , = d»A v - d v A» , (3.4) 

and a strangeness-conserving term of O(e 2 p ) [ p^j . 

e 2 F A Z(Q L Q R ) . (3.5) 

The numerical value of the parameter Z can be determined from the mass difference of 
charged and neutral pions. The relation M 2 ± — M 2 = 2e 2 ZF 2 implies Z ~ 0.8. 

Finally, we have to introduce a weak-electromagnetic term characterized by a coupling 
constant g ew k, 

e 2 F A (TQ R ) , T = g cwk G 8 F 2 A + h.c. . (3.6) 

Note that to lowest order only a single (linear independent) term of this type can be 
constructed once the relation 

Q L A = AQ L = ~A (3.7) 
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is taken into account. This term is the lowest-order chiral realization of electromagnetic 
penguins || [l7j and transforms as (8^, 8_r) under G when the Q R spurion field is "frozen" 
to the constant value ([2.7]) . By chiral dimensional analysis we expect the coupling constant 
fl'ewk to be of 0(1). A recent estimate in Ref. |2| corresponds in fact to g cw k = —1.0 ± 0.3 
(see also Ref. @). 

Summing up all these contributions, our lowest-order effective Lagrangian assumes 
the form 

A = ^ (u^ + x+) + F 2 (EuX) 



-F^ + e 2 F 4 Z(Q L Q R ) + e 2 F 4 (TQ R ) . 



(3.8) 



Using (|2.13|), one easily verifies that (|3.8|) is CPS invariant. 



4 One-loop divergences 

For the construction of the one-loop functional, we first add a gauge-breaking term (we 
are using the Feynman gauge) and external sources to ( |3.8| ): 

c 2 ^c 2 - l -(d^f-j^ . (4.i) 



Then we expand the lowest-order action associated with fl4.1|) around the solutions ip c \, A^ 
of the classical equations of motion. In the standard "gauge" u R (ip c {) = ul(<P c \V ='■ u(ip c i), 
a convenient choice of the pseudoscalar fluctuation variables & (i — 1, ... ,8) is given by 

u R = u cX e^ 2F , u L = ule-^ 2F , fcfad) = , (4.2) 

with the Gell-Mann matrices \ (i — 1, . . . , 8). The photon field is decomposed as 

A» = A£ + e" (4.3) 

with a fluctuation field £ M . In the following formulas, we shall drop the subscript "cl" for 
simplicity. The classical equations of motion take the form 

V/ = t -{ X --l(x~))+2te 2 F 2 Z[Q R ,Q L ] 

+i[«X,S] - 2(V„K,H} - |<V M K,S}» 
+2 i e 2 F 2 [Q JJ ,T] , (4.4) 

eF 2 

□A M = J fl + —-{u tl {Q R -Q L ))+eF 2 (E{Q R -Q L ,u^}) , (4.5) 



where 



v M = <9 M + [r M , 
i 

2 



*V = -[^{d^-irju + uid^-il^) . (4.6) 



The solutions of (fOj) and (|4.5|) are uniquely determined functionals of the external 
sources t> M , a^, x^ J^- (Note that the usual Feynman boundary conditions are always 
implicitly understood.) 

Expanding (|4.1|) up to terms quadratic in the fields £ i5 e M , we obtain the second-order 
fluctuation Lagrangian £( 2 \ The one-loop functional Wl=x is then given by the Gaussian 
functional integral 

e w ^ = J[d^de ll ]e i S^ . (4.7) 

In our case, reads 

£ (2) = ^(v^v^ + ^«-^(«X + x+)e 2 ) 

+ JF*Z(ZQi£Q R -±e{Q L ,Q R }) 

+ ^mv»^ - !{«„, K,e 2 }} + - 2i[e,s]{v^,^}) 

+ e 2 F 4 (T(eQ^-^{^ 2 ,Qi ? })) 

+ ^de" + ^((Qr ~ Ql?)s^ + e 2 F 2 (Z(Q R - Q L f)e^ 

ie F 2 p F 2 

(K, Q R + Q L ]0^ + V((2« - QzW^e* 



4 u ^ " ^' 2 
ieF 2 ,_ rr ^ ^ ^ u zeF 2 



^ ■(S{[Q /? + Q L ,e],M)^-^-([^S]{Q i? -Q L ,^})£ 
+ eF 2 (S{Q K -Q L) V^})^, (4.8) 

where 

f = . (4.9) 

In the next step, we perform the field transformation 

£-f-{E,£} + !<S£>l. (4.10) 

Because of (A) = 0, we do not pick up an additional contribution from the Jacobi 
determinant and the fluctuation Lagrangian (|4.8| ) assumes the form 



where 

<■ • • (4-12) 

The explicit expressions for 7^, cr^-, k, af , 6j are given in the Appendix. 
The divergent part of the one-loop functional, 

Wfl x = J d d xCt =1 , (4.13) 



is determined by 

££i = - , A J d _ 4) M^-r + ^ 2 ) - <«* + <W)< 

+ o( 6 < 6 *) 2 - - ^ + 2ft;2 ] > ( 4 - 14 ) 

where 

7^ = - <%7 M + [7m, >] • ( 4 - 15 ) 

This formula can easily be derived from the well-known second Seeley-deWitt coefficient 
for bosonic systems |19| . 



5 The chiral Lagrangian at next-to-leading order 

We are now in the position to construct the most general local action at next-to-leading 
order which will also renormalize the one-loop divergences discussed in the previous 
section. 

The strong part of the local action of 0(p 4 ) is, of course, generated by the well-known 
Gasser-Leutwyler Lagrangian Jl3| associated with the low -energy constants L\, . . . , L\2 ■ 
In the presence of virtual photons, the structure of the operators given in |13j remains 
unchanged. The only necessary modification is the inclusion of the dynamical photon field 
in the generalized "sources" and r M (see ( |2.3| )). The divergences corresponding to the 
strong sector of ( |4.14| ) are absorbed by the divergent parts of the Li [13]. In the relevant 



case of chiral 577(3), the strong terms generated by ( 4.14|) can be written immediately as a 
linear combination of the 0(p A ) operators of the Gasser-Leutwyler basis without using the 
equations of motion ( |4.4[) or ( |4.5|) . Consequently, no additional (weak-) electromagnetic 
terms are induced at this point. 

The strangeness-conserving terms of 0(e 2 p 2 ) have been constructed by Urech ||. His 
list of electromagnetic counterterms is associated with the coupling constants K\, . . . , K14. 
In this case, ( |4.14j ) leads to that canonical basis only after the use of the equation of motion 
( O ). In this way, also some divergent weak-electromagnetic contributions of 0(G$e 2 p 2 ) 



are generated. 

For the octet part of the nonleptonic weak Lagrangian of 0(GfP 4 ) [20] we refer to 



the standard form of Ecker, Kambor and Wyler || with couplings N\, . . . , ./V37. Again, 
because of the mismatch between ( 4.14Q and the standard basis, the equation of motion 



has to be used and the (purely) electromagnetic piece in ( |4.4j ) induces divergent terms of 
0(G 8 e 2 p 2 ). 

Finally, we have to construct the most general weak-electromagnetic Lagrangian of 
0(Gse 2 p 2 ). Some parts of this Lagrangian have appeared before in the literature |], |], ^TJ . 
The complete minimal Lagrangian of 0(Gse 2 p 2 ) takes the form 

32 

£ G8 eV = G 8 e 2 F 4 ]T ZiQ t + h.c. , (5.1) 
i=i 



7 



with operators Qi of 0(p 2 ) and dimensionless coupling constants Z$. A linear independent 
set of operators is given by 

Q 1 = (A{Q R , X+ }}, 

Q2 = (AQ R )(x+) , 

Q 3 = (AQ R )( X+ Q R ) , 

Q 4 = (A x+ )(QlQr) , 

Q 5 = (AuX) > 

Q 6 = (A{Q R ,uX}) , 

Q 7 = {AuX)(QlQr) , 

Q s = (AuA(QX) , 

Q 9 = {Au»){Q R u») , 

Qio = (A^)({Q l ,QhK) , 

Qn = (A{Q R ,^})(Q L ^) , 

Q12 = (A{Q fl ,« M })(Q fl «") , 

Q 13 = (AQ R ){uX) , 

Qu = (AQr){uXQr) , 

Q 15 = (AQ fl )(«X(Qi-Qfl)) > 

Qie = (A X+ ) , 

Q17 = |(Ax + > + (A{Qfl,x+}> + <A[Qfl,X-]> , 

Q 18 = {A{Q r , X+ })- 1 -{A[Q r ,X-]) + {HX-QlQr-QrQlX-)) 
- ^(AQ R )( X+ ) - (AQ r }( x+ Qr} + (A x+ }(QlQr) , 



Ql9 


= (AQ r )( x+ (Ql-Qr)) , 


Q20 


= i((V M A)[Q L ,^]> , 


Q21 


= i((V M A)[e Jl ,«'*]> , 


Q22 


= i((V„A)(QXQ fl - , 


Q23 




Q24 


= i(A{u^Q L )Q R -Q R {^Q L )u,)) 


Q25 


= i{(V^A)(u"Q R Q L - QlQru*)) , 


Q26 


= i(A{u,Q R {V»Q R )-{V»Q R )Q R u,)) 


Q27 


= t(A(Q R u^Q R ) - (V^Qr^Qr)) 


Q28 


= <(V„A)(V"Q L )> , 


Q29 


= ((V^KVCh)} , 


Q30 


= (A(V»Q R )(V»Q R )} , 


Q31 


= ((V m A){V m Ql, Qr}) , 
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Q 32 = ((VAHV^Q*, Q L }) , (5.2) 

where 



V^A = V„A + Uup, A] = u(D„X)v^ 

i 



V m Qk = V m Qk-|[^,Q k ] =M t ( J D M g i? ) M , (5.3) 



with 



D^Ql = 9nQ L — i[lfi, Ql) , 

D^Q R = d„Q R -i[r„Q R \ . (5.4) 

For the construction of the list of local terms (|5.2|) we have used CPS invariance, the 
relations (|3.7]) and 

Ql,n = l± + IQl,r , (5.5) 

the Cayley-Hamilton formula, partial integration and the equations of motion (|4.4j) . 

If the spurion fields Ql,r and A are fixed to the constant values in ( |2.7| ) and ( |2.11| ), 
respectively, then C Gae 2 p 2 transforms under G as 

(8 L , 1 b ) + (8 L , B r ) + (27 L , 1 R ) + (27 L , 8 R ) + (8 L , 27 R ) . (5.6) 

This structure is richer than the one of the 0(G&) terms in C 2 and also of the weak 
four-fermion effective Hamiltonian [fT2]| . The last two pieces, in particular, which are 
responsible for AI = 5/2 transitions, have no analog in the effective Hamiltonian of 
dimension six. 

The operator Qi6 does not contribute to on-shell matrix elements [ID, ^2], |23|, 20] . The 



terms Qn, Qi$, Qig vanish for electrically neutral (pseudo) scalar sources, 

[X,Q} = 0, (5.7) 

which is, of course, the case for all realistic physical processes. Also the operators 
Q20, ■ ■ ■ Q32 are irrelevant for practical purposes. Because of ( |5.3[ ) and Q5.4p, they con- 
tribute only in the presence of non-vanishing external (axial-)vector sources. 

The coupling constants Z±, . . . , Z 12 appear in the amplitudes of K — > 2n decays. 
The operators Q13 and Q14 do not contribute to K — > 2tt but they enter for K — ► 3n. 
Qi 5 involves at least five pseudoscalar fields and is therefore irrelevant for K decays. 



A few linear combinations of the operators in ( |5.2j ) were already given some time ago 
by de Rafael 0. His list was restricted to terms contributing to K — > 2ti, neglecting 
contributions ~ M% and those renormalizing Gg. A more recent extension of de Rafael's 
list can be found in Ref. ||. However, their Lagrangian is still incomplete even for the 
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K — > 2n amplitudes, as we shall discuss in the following section. There is in addition an 
obvious misprint in the operator multiplied by Sq in |J, which would be in conflict with 
chiral symmetry. Some of the operators in ( |5.2| ) have also appeared in attempts |H[ to 
bosonize the AS = 1 four-fermion effective Hamiltonian. 

The low-energy couplings Zj are in general divergent. They absorb the divergences of 
the one-loop graphs via the renormalization 

Zi = ZT(ii) + ZiA(ii) , 2 = 1,..., 32, 
AM = |^{^4-^H4vr)+r'(l) + l]}, (5.8) 

in the dimensional regularization scheme. The coefficients zi, . . . Z32 are determined in 
such a way that the divergences generated by ( (4.14 ) are cancelled: 



3 7Z 
1 

~2 ' 
3 

2 ' 

9 3 
1 

"2 ' 

0, 
1 

"2 ' 

§ + 6Z. 

(5.9) 

As already discussed above, the values in this list depend on our conventions for the 
basis systems in the strong, electromagnetic and weak parts of the next-to-leading order 
Lagrangian. The Zi given in (|5.9|) have to be used together with the divergent parts of the 
coupling constants Lj [jnj, K{ and iVj M, respectively. The divergences involving the 
electroweak penguin coupling g ewk are independent of this choice of basis and they agree 



Zl = 


17 
~12 ~ 


3Z + — (?cwk > 


22 = 


1 u_ 16 7j _ 

- 1 + — ^ + S'ewk , 


z 3 


-\ + iz, 


24 = 


Z5 = 


-2 , 




^6 = 


= - + 5Z + -g ew k , 


z 7 -- 


=| + 5Z, 


^8 = 


z 9 = 


11 

nr + 


-Z + 2p C wk , 


zio 


2 


Zll 


= -?-2Z, 

2 


Zl2 = 


z 13 = 


35 
~12 ~ 


- 3Z + (? ew k , 




= 3 + 15Z , 


z 15 


= ^ + 15Z, 


z ie — 


zn = 


1 + 2Z, 


z 18 


= \ + 3Z, 


z 19 


= 4Z, 


z 20 — 


Z21 = 


1 

6 ' 




Z22 


= 3 + 6Z , 


Z23 


= -3-9Z , 


Z 2A = 


Z25 — 


-3Z , 




z 26 


= -1 , 


%27 


= , 


Z 2S — 


Z29 


1 

"2 ' 




z 30 


= 0, 


Z31 




z 32 — 



with a recent calculation of Cirigliano and Golowich | 23fl . Note that g ew k appears only in 
the couplings of (8^,8^) operators. This is because the lowest-order term proportional 
to <?ewk is already of 0(G$e 2 ). Therefore, the 0(Gse 2 p 2 ) terms proportional to g ew k arise 
from the product of the lowest-order (8l,8r) weak operator times the 0(p 2 ) invariant 
part of the strong Lagrangian. 

The renormalized low-energy constants Z\{p) are in general scale dependent. The 
coefficients Zi govern this scale dependence through the renormalization group equations 

= ■ (5 ' 10) 
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By construction, the complete generating functional at next-to- leading order is then scale 
independent. 



6 K 



7T7T 



In the modern framework of chiral perturbation theory, electromagnetic corrections for 
K — > txtx decays to 0(G 8 e 2 p 2 ) were discussed by de Rafael || and have been treated in 
more detail by Cirigliano, Donoghue and Golowich [§, H ||. Together with corrections of 
O (G$(m u — md)p 2 ) ||], the complete isospin-breaking effects of next-to-leading order have 
obvious phenomenological implications, from the AI = 1/2 rule to CP violation [ [25| . 

In this section, we present the tree-level contributions to the K — > txtx amplitudes from 
the Lagrangian (|5.1|) . We compare those amplitudes and in particular their divergent 
parts with the results of Ref. ||. Using our own one- loop calculation of isospin-breaking 
corrections [0 and the heat-kernel results (|5.9|) , we find that the complete amplitudes of 
0(Gse 2 p 2 ) are indeed finite. We demonstrate the cancellation of divergences explicitly 
for the subset of amplitudes proportional to the electromagnetic penguin coupling g ewk 
defined in ( |3.6| ). 

From the Lagrangian ( |5.1| ) of 0(G$e 2 p 2 ), we obtain the following amplitudes in units 

of C owk := iG 8 e 2 F: 

■ C cwk V2[(M 2 K -M 2 )(2Z 1 + AZ 2 

-2/3Z 7 ) + M 2 (6Z X + 6Z 2 - Z 6 ) 

= C ewk V2(M 2 K -M 2 )(-Z 5 + 2/3Z 6 
-2/3Z u - 2/3Z 12 ) , 



ewk 




TX TX 



4/3Z 3 + 4Z 4 - Z 5 - 1/3Z 6 



A(K° 
A(K + - 



7X°7X°) 



2/3Z 7 + Z 8 + Z 9 + 2/3Z 
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7T 7T 



a 



ewk 



{M 2 K - M 2 )(2Zx + AZ 2 - 4/3Z 3 - Z 6 - Z s - Z 



-4/3Z n - 4/3Zi 2 ) + M2(6Zx + 6Z 2 - Z 6 



2/3Z 10 
(6.1) 

0. In addition, the 



These amplitudes agree with Ref. [g] for Z 3 = 3Z 2 , Z 10 = Z\\ 
coefficients s 8 , s 9 in Eq. (35) of should be multiplied by 2/3. 

In the SU(3) limit for the mass matrix (|3.2j) , the amplitudes ( |6.1| ) satisfy the relations 



A(K° - 
A(K° 



IX TX 
> 7T 7T 



)SU{3) 
)SU(3) 



V2A(K+ ^ 7X + 7X°) SU{3) , 

, (6.2) 
7X7X transitions in the presence of electro- 



in accordance with a general theorem on K 
magnetism [0. 

The divergent parts of the Zi in Q5.9D give rise to the following divergent tree-level 
amplitudes, with A(/x) and Z defined in ( |5.8| ) and (|3.5|), respectively: 

CewkV^A(^) [Mi (-3 - 27Z + 13/2^ ewk ) 

+M 2(-3 + 36Z + 7(7 cwk y , 

C ewk V2A(fi)(M 2 K - M 2 )(2Z + 3 9cwk ) , (6.3) 



A(K° -> 7T+7T-) 



div 



A(iT° 
A(K+ - 



JJ\ 



C ewk A(/i) [M|(3Z + 7/2(? cwk ) + M 2 (-Q + QZ+ 10g cwk ) 
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Figure 1: Loop diagrams for K — > -kit involving the electromagnetic penguin coupling 
g ew k- The associated vertex from the Lagrangian ( |3.6| ) is denoted by a crossed circle. 
Normal vertices are from the lowest-order strong Lagrangian (|3.1| ). 



The (ultraviolet) divergences in ( |6.3|) arise from three different sources: 

• Photon loops proportional to G$e 2 ; 

• Loops involving the electromagnetic coupling (|3.5| ) proportional to Gge 2 Z; 

• Loops involving the coupling fl3.6|) proportional to Gge 2 g ewk . 



Strong and electromagnetic wave function renormalization [3f| is included in all three 
categories. 

We have performed a complete calculation of K — > tttt amplitudes to 0(Gge 2 p 2 ) and 
0(G 8 (m u — rridjp 2 ) 0. For m u = m^, we find that the explicit loop divergences are exactly 
cancelled by the divergent tree- level amplitudes (|S~3|) . We exhibit those cancellations in 
detail for the divergences proportional to g cvlk . Divergences arise both in loops with 
an electromagnetic penguin vertex shown in Fig. [I] and from (strong) wave function 
renormalization of tree diagrams from the Lagrangian fl3.6| ). 

In the exponential parametrization, the divergences due to the diagrams of Fig. |l]take 
the form 

^/2 

A(K° -> 7r + 7r-) loops = -^-C ewk g cwk A(^)(7M 2 K + 8M 2 ) , 

A (K° ^ 7rV) loops = -3V2C cwk9cwk A(fi)(M 2 K - M 2 ) , 
A (K + -> 7r + vr°) loops = -C cwk9cwk A(fi)(M 2 K /2 + 7M 2 ) . (6.4) 

Wave function renormalization (again in exponential parametrization) leads to 

A(K° -> 7T+7r-) wfr = -3V2C cwk g ewk A(fx)(M 2 K + M 2 ) , 
A(K° - ttVU = 0, 

A (K+ -> tt+ttV = -3C ewk ^ ewk A( M )(M| r + M^ . (6.5) 
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The sum of ( p.4| ) and (|6.5|) is parametrization independent and it is exactly cancelled by 
the terms in (|6.3| ) proportional to g ew k- 

We have exhibited (part of) the loop divergences explicitly also because we do not 
completely agree with the results of Ref. [|J. Although the divergences due to photon loops 
are identical, we obtain different results for some of the other divergences^]. Only for the 
channel K° — > 7r + 7r~, there is complete agreement for all three types of divergences. 

The complete amplitudes of 0(Gse 2 p 2 ) and 0(G$(m u — md)p 2 ) together with a phe- 
nomenological analysis will be presented elsewhere |7j. 



7 Conclusions 

We have supplied the missing ingredients for a complete analysis at next-to-leading order 
of the combined strong, nonleptonic weak and electromagnetic interactions of mesons. 
The main results are: 



The complete and minimal Lagrangian (5.1) of 0(Gge p ) contains 32 operators Q 



and associated dimensionless coupling constants Zj. Of these 32 operators, only 14 
are of immediate phenomenological relevance. We have ordered the terms in a way 
most suitable for applications: the first 12 operators contribute to K — > 27r decays 
whereas the remaining two enter in K — > 3tt amplitudes. 

ii. The one-loop divergence functional ( |4.13|) determines the renormalization of the 
effective theory. Together with the previously known divergences, the new terms 
( |5.9| ) in the coupling constants Zi ensure that the complete amplitudes for strong, 
nonleptonic weak and electromagnetic interactions of mesons at next-to-leading 
order are finite. 

As a first application, we have presented the tree-level amplitudes of 0(G$e 2 p 2 ) for 
K — > 7T7T decays. The associated divergent parts cancel with the explicit one-loop diver- 
gences to yield finite and scale independent decay amplitudes. 
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1 V. Cirigliano has informed us that they now agree with the divergences ( |6.3[ ); see forthcoming erratum 
for Ref. @. 
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Appendix 



The quantities occurring in (|4.14j) can be decomposed with respect to (explicit 2 ) powers 
of e and Gg in the following way: 



<Tij 

7„ 

K 



(7 ij\e°G° + Cr ij\e 2 G° 8 + a ij\e°G s + cr «ile 2 G 8 

AM i AM 

a i leGg a i |eG 8 > 

^leG° + h\eGg, , 

K \e 2 G° + K |e 2 G 8 • 



(A.l) 



The explicit expressions for the various terms are given by 



a ij\eOG% 
a ij\e 2 G% 
a ij\e°G 8 



a ij\e 2 Gg 



TylePGS 



7ij|e°G 8 



AM 
°i leGg 



g{(«/ + X+){A i; A,}) - -(UfiXiU^Xj) , 
e 2 F 2 Z(i{Q /? , Q L }{A,, A,} - XiQ R XjQ L - X 3 Q R \Q L ) , 
-({H, AiKAX) + ^({S, AyKA^) 

- ~<(«X + X + )({A l; {5, A,}} + {A„ {5, A,}})) 
+ -(EX^ix+Xj) + -(EXj)(x+Xi) 

+ J<SK,K,{A i) A i }}}> 

- ~<H(K, A 4 }K, A,} + K, X 3 }{u», A,})) 

+ ^(K, VS]{A,, A,}) + ^<[V%, ~]{A l; A,}) 

- ^{{Ai, A^VS) , 

e 2 F 2 Z (E(XiQ,RXjQL + \ s Qr\Ql + Qr^QlXj + Q R \jQ L K 
+ X t Q L XjQ R + XjQ L XiQ R + Q L X t Q R X 3 + Q L XjQ R X t 

- A,{Ql, Q r }X 3 - Xj{Q L , Q R }X t - ^{{Q L , Q R }, {X t , A,}})) 
+ e 2 F 2 (T(i{{A,, A,}, Q R } - X t Q R X 3 - X 3 Q R Xi)) , 
-^HA,,A,]> , 

ia^sjK.A^-i^.siK.A,}), 



(A.2) 
(A.3) 



(A.4) 



ieF 



u»[Q R + Q L ,\i]) 



(A.5) 
(A.6) 
(A.7) 
(A.8) 



'Note that e also appears in the vielbein (2.2) and in the connection T M (4.6) via (|2.3|). 
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a i \eGg 



bi\eG° 8 
h\eGg 



K \e 2 G° 8 



7Pr P W 

— (-(u^KQn - QrX^)) + —{E^KQl - Q L \iU»)) 

ip W ?p F 

- —(E(u»Q R \i - KQr^)) + T (H(A,/Q fi - Q R u»\)) 

ieF pF 
+ T (H(Wfii - Ql^K)) - —{(Qr - Ql){V^H, A,}) , (A.9) 



^■{{Qr-Ql)\) , 

e -^(E{Q R -Q L ,\}) 
e 2 F 2 



e 2 Gs 



■((Qr-Ql) 2 } , 
2e 2 F 2 (E(Q R — Ql) 2 ) 



2 

.2 772, 



(A.10) 
(A.H) 

(A.12) 
(A.13) 



The expressions ( |A.12 ) and ( [A. 13|) are included for completeness only; k does not 
contribute to the order we are concerned with. 
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